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Abstract 

We consider the QCD sum rules approach for A-isobar magnetic 
form factor in the infrared region < Q 2 < lGeV 2 . The QCD sum 
rules in external variable field are used. The obtained form factor is 
in agreement with quark model predictions for the A-isobar magnetic 
moment. 



1 Introduction 



The QCD sum rule method suggested by Shifman, Vainshtein and Zakharov 
(SVZ) in the pioneering paper [|lj] becomes now a universal tool for calculating 
different properties of low-lying hadronic states. Using the original version of 
this method, the meson Jl| and baryon [0 masses were found from the sum 
rules for two-point correlation functions. Using the three-point correlation 
functions, hadron form factors at intermediate Q 2 can be obtained ||. Un- 
fortunately, this method does not work if one tries to calculate form factors 
in the infrared region < Q 2 < lGeV 2 due to power corrections 1/Q 2n at 
Q = 0. The new method - QCD sum rules in external field was suggested 
in ||, and using this method nucleon magnetic moments were found as 
well as baryon axial couplings ||. Then this method was formulated for a 
variable external field [0] which gives a possibility to calculate form factors 
at QV 0. 

In [0 we have formulated a new method for calculating hadronic form 
factors in the infrared region. To study a form factor at nonzero Q 2 , it is 
necessary to introduce a variable external field. The calculation of a polar- 
ization operator in this field encounters a number of difficulties as compared 
with the case of a constant external field. The arising problems and methods 
to avoid them were discussed in detail in the paper 0. 

Let us note that in the papers || and |§, the pion and nucleon charge 
radii were considered, using the methods similar to ours. However, the results 
obtained in || are connected not with the calculation of the total form factor 
F(Q 2 ), but only the first derivative at zero momentum transfer < r 2 >~ 
F'(Q 2 )\q =0 . In || pion form factor was considered. The aim of the paper is 
to use the general method for calculating A-isobar magnetic form factor in 
the infrared region. 

2 Polarization Operator in Variable External 
Field 

To compute A-isobar magnetic form factor we shall consider the following 
correlator in an external variable electromagnetic field: 

n£(p,*0 = i J 27c5(kx)e ipx d 4 x < T{ Vti (x),fj u (0)} > v (1) 
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+oo 



where 



dzi J e i{p+kz)x d 4 x < r{^(ar),^(0)} >y 

/+oo 
dzriL(p + kz,k) 
-oo 

n^(p, fc) = i / e^x < T{^(x), 77,(0)} > y (2) 

^(a;) = e abc (u a C^u b )u c (3) 

is the quark current with the A-isobar quantum numbers suggested in the 
first paper in Ref.@, u is the u-quark operator; a , b and c are the color 
indeces, e abc is the antisymmetric tensor and C = —C T is the charge conju- 
gation matrix; index V means the vacuum average in the presence of weak 
external electromagnetic field that is responsible for adding to Lagrangian of 
the following term 

AC = -V^ kx J2e f qf(x) lfl q f (x) = -V^tfe** (4) 
/ 

where ej is the charge of the quark with flavor /, and are the ampli- 
tude and the momentum of the classical external field. This correlator was 
suggested in @. 

Now let us discuss the reason why we need to introduce the ^-function in 
the correlator (p. To calculate the correlator @ at p 2 ~ — lGeV 2 , (p + q) 2 ~ 
— lGeV 2 , k 2 = —Q 2 < 0) we use operator product expansion in the presence 
of external variable field So we need to know nonperturbative quark 
propagator in the field 

<:T{q^x),q b (0)}:> v (5) 

=< T{q a a (x),q*(0)} > v - < T{q" a (x),q*(Q)} 

Where :: denotes a subtraction of perturbative contribution. It is possible 
to find perturbative part of this propagator in the form of expansion over the 
coupling constant. To take into account nonperturbative interaction of the 
quark with the external field, we expand eq.(Bp over x^ 

<: T{q a a (x), q b (0)} :> v =<: q«(0), g$(0) :> v (6) 
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+x p <: £>^(0),^(0) :> v + -x p x u <: D„D„g£(0), q b p (0) :> v +... 

It is clear that the n'th term of expansion (§) can give dimensionless factor 
(kx) n . Effectively it means that highest terms contribution of expansion (^) 
into a polarization operator (fj) is not suppressed because this factor {kx) n 
corresponds to the factor (kp)/p 2 ~ 1. To kill the dangerous contributions 
of the terms ~ (kx) n we insert 8{kx) into the correlator ([I]). This correlator 
([!]) can be calculated in a form of series over 1/p 2 . Therefore at respectively 
large — p 2 this correlator can be calculated with a good accuracy using only 
the first few terms in expansion 

The nonperturbative quark propagator in the external field has the fol- 
lowing form: 

5 ab ~ 

<: <(x),4(0) :>= e u — {V^Ii^k 2 ) 

+ ^ pX )^k p V x U 2 (k 2 ) + (Vx)[i <: W> :> +h 2 U 2 (k 2 )}5 a p (7) 

+ l -V aP {kx)k 2 Tl l {k 2 ) + l -{a pX ) a pk p V x {kx)Tl 2 {k 2 ) 

~x p e pup x{^\^) a pk u V p k 2 T[i(k 2 ) 

+-^x 2 (a pX ) af3 k p V x {<: W :> -ill? (k 2 ) - 2U^(k 2 ) + ik 2 Tl A {k 2 )} 

+ (Vx)(kx)k 2 U 3 (k 2 )5 a(3 + ^(a pX Uk p V x U 4 (k 2 ) 

<= W ->o +nf (A; 2 ) + 3^ 2 n 3 (£; 2 ) + ill^k 2 ) 

-^ 2 n 4 (A; 2 )] + ^- Xp (a pu )^k u [^ <: ^ :>o -Ilf (A; 2 ) 

-iR^k 2 ) + 3^ 2 n 3 (A; 2 ) - yn 4 (A; 2 )] 
+ {terms with an odd number of ^ — matrices)} + 0(x 3 ) 
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The correlators Ui(k 2 ) are defined as follows: 

ejli(k 2 )(k 2 g py - k p k u ) = 

i f e ikx d 4 x <: T{^q f 'y lt q f (x),u'y v u(0)} :> 
/ 

e u Il 1 (A; 2 )(A;^ p - k v g pp ) = (8) 

i I e ikx d 4 x <: Ti^QflpQfi^^a^uiO)} :> 
/ 

e u Il 3 (k 2 )(k 2 (g pp k u + g vp k p ) - 2k p k v k p ) = 
i I e ifcx 'rf 4 x <: T{52e f q f 7 P qf(x),u D {il D v} u(0)} :> 

c^n^A; )k p k L ,(g pr k p gppk r ) + ... = 
f e ifcx 'd 4 x <: e f q f j p q f (x), w pr D^D^y w(0)} :> 

(fc )(k p g pu k u g pp ) = 
2 / e^ 4 x <: T{^e / g /7p g / (x),^G'^r U (0)} :> 

){kp9pK kvgpp) = 

e lk *d 4 x <: r{Ee/?/7pff/(a:),ftfiG^*"7fett(0)} :>o 
/ 

Perturbative contributions are subtracted in correlators (|8|). These expres- 
sions were obtained in 0. 

In this paper we neglect operators U G (k 2 ) because their contribution into 
a sum rule is small (see 0). 

3 The Sum Rules 

In this Section, we obtain a sum rule for the A-isobar magnetic form factor. 
First, we should choose a tensor structure which has a contribution from the 
magnetic transition between baryon states with quantum numbers J = 3/2. 
To this end, we consider the contribution of two baryons with masses mi and 
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rri2 into the polarization operator Tl v (p, k) (|2|) 

V p < | 7fc | Ai >< Ai | 3 e P m I A 2 >< A 2 | ^ | > 
(p 2 — m±)((p + A;) 2 — m 2 ) 



(9) 



where Ai and A2 are baryon states with masses mi and 777.2 respectively. Here 
we consider the case when only spinor parts of the Rarita-Shwinger fields in- 
teract with a photon. In such case, the matrix element of the electromagnetic 
current has the following form: 

< N, I f p m I N 2 >= 

v^(p)g,AUk 2 )i P + igS^A + iM* 2 )*,]*^ + k) = 

4\p)9AUi2{k 2 ) + Vi2(k')) lp + V P ^Q + ^mk p ]v^\p + k) (10) 
<0\ Vfl \N, J p = 3/2+ >= Xv p (p) 
VfM = Pti + (p + 

where v^{p) is a Rarita-Shwinger spin- vector satisfying the Dirac equation: 
(p - m> M (p) = 0, 7^ = 0, p^ = 0. 

Using ( prp| ) we can transform (|^) to the following form: 

A i A2 u l a a n fi* 



V,[g^pa p p 2 G M (k 2 )/3 (11) 



(p 2 — 777 2 )((p + k) 2 — ml) 

+ (other structures with 7 M placed at the beginning and 7^ at the end of them)} 

where Gm is the magnetic form factor. 

It is important to note that there is no spin-1/2 baryon contribution in 
the structure gp, v P\lpP2 which has the following amplitude: 

< I n p I J = 1/2 >= (Ap M + B lp )u(p) 

where (p — m)u(p) = and Am + 4B = 0. 

From (0), it is obvious that the structure g^pij^ (where pi = p and 
P2 = p+k) contains magnetic transition only, since G^/3 = /i2(^ 2 )+V 9 i2(^ 2 ), 
where G 1 ^ is the magnetic form factor. So, we shall further consider only the 
structure g^PxlpVi- 

Now let us discuss the factor 1/3 which have appeared in eq . (|TT|) . Con- 
sider interaction of spin-3/2 particle with the electromagnetic field: 

a(* M (P 1 )^(P 1 + P 2 ) P * U (P 2 ) + ib^ p {a pX g pu /2 + 2g lv g vX )^ v F pX (12) 
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where \o p \ = [jpjx], ^ is Rarita-Shwinger spin- vector field ( (P — m) — 
0,7^^1/^ = 0), F p x = d p A\ — d\A p . The first term of eq. (|T2|) corresponds to 
the spin-independent part of electromagnetic interaction of |-spin particle 
and the second one describes the spin-dependent interaction. To express the 
value of the magnetic moment (at Q 2 = 0) let us consider the case when 
A Q = 0, P = m,F p x = dpiSxjeijkHk, where H k is magnetic field. Then we 
have 

m m {g mn&ijtijk 

)V n H k (13) 

where S fc = diag(a k , a k ) Now we see that the operator (£ fc g mri + 2ie mnk ) is 
equal to 2S k where S k is spin operator for the Rarita-Shwinger field. So, the 
maximal energy of the particle in the magnetic field is equal to 

E int = 3bH = liH, (14) 

where /i, by definition, is magnetic moment or magnetic form factor at Q 2 = 
0. Thus, we have 

/i = 36 (15) 

Now let us consider the double dispersional relation for the function at 
tensor structure g p V V\l pV'iY p- 

f{ pi pi q 2 ) = Lrr '|'--?ff> . , del 

n 2 Jo Jo {Pi + si)(Pr f + s 2 ) 

where P 2 = —p\, P 2 2 = — p 2 , Q 2 = —k 2 > and p(si,s 2 ,Q 2 ) is the spectral 
density. Due to reasons mentioned above, we cannot calculate p(pi,P2,Q 2 ) 
directly, but we can consider the double Borel transformed structure function 
of correlator ([]]). Notice, that under replacement p\ — > p\ + kz, p 2 — > p 2 + kz 
the structure pi^p^V^ transforms to 

PilpfaVp -> (pi + kz)V{p 2 + kz) = piVp 2 - 2{Vpi)piz - (17) 
2{Vp 2 )p 2 z + {p\ + p\)zV + z 2 {2{Vk)k - Vk 2 ) 



and all other structures in ( |TT|) with smaller number of 7— matrices near 
g^ u , could not be transformed into p\Vp 2 . So we can extract gp V piVp 2 
in the integral ([!]). Then the integral representation for the double Borel 
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transformed structure function under consideration is obtained from (|TJj) by 
applying the operator O: 



r+oo 

Of {Pi P 2 2 ) = / dzBftBptfdP! + kzf, (P 2 + kzf) (18) 

where 



n! V 5p2 / 

P 2 /n = M 2 

M 2 is the Borel parameter. 

Applying the operator O to the left-hand and right-hand sides of double 
dispersion relation ( |16D we get 

1 r +l/2 n2 z 2 

roo roo Sl + S2 ( S1 -S 2 )Z 2 

Jo Jo dsids2ex v(~^p- + )P( S 1^2,Q) (19) 

where M 2 = M 2 = 2M 2 . 

We shall use the sum rules (|19|) to calculate the nucleon form factor. The 
function f(M 2 ,Q 2 ) is calculated using operator expansion in the external 
field, and the spectral density p(si, S2, Q 2 ) is saturated by intermediate states 
with quantum numbers of the current (|3|). 

There are two different types of intermediate state contributions into ([L9|) . 
The first is responsible for the diagonal transitions between the states with 
equal masses. The second is responsible for nondiagonal transitions between 
the states with different masses. In the first case the right-hand side of (|T^) 
obviously has the form 



A 2 ^(Q 2 )e- mVM2 r l2 e^dz (20) 
3 .7-1/2 



where A 2 is the square of the residue of the state with mass m into the 
current rj^ defined by formula (|T0|) , Gm(Q 2 ) is the corresponding magnetic 
form factor. 
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In the second case for the transition between states with masses mi and 
m 2 we get: 

„ m 2 +m 2 /.+1/2 Q 2 z 2 + (m 2 -m 2 ) z 

A 1 A 9 ^(Q 2 ) e -^7^ / e s^^dz (21) 

3 7-1/2 

Now, to investigate the properties of ( BPp and ( |2"T| ) let us expand them in the 
series on Q 2 /M 2 

2 G M {Q 2 ) m 2 /M 2 1 Q 2 1 g 2 2 
A 3 li+ 12M 2 + 5-2^M 2j + - j 



? 2M 2 , ,m? - mi 
75 =■ smhl — 

m? - m| V 2M 2 



A^gjg(Q 2 )e-^ „r"„ 2 sinh( ""o^" a ) >< (23) 



Q 2 1 M 2 m 2 -m 2 2M 4 

+ M 2 4 _ m 2 -m 2 ( 2M 2 } + (m 2 -m 2 ) 2j + 



From (|22| ) and (|23f) we see that diagonal transitions of the excited states are 
exponentially suppressed compared to the A-isobar contribution in Q22]). Let 
us write the non-suppressed part of the contribution from the nondiagonal 
transition between the nucleon with mass and an excited state with mass 

AAA* 2 21 1 + 

3 m A » — m A 
Q 2 1 M 2 m^-m 2 2M 4 1 

M 2 4 mi, - ml [ 2M 2 } + (m 2 A , - m 2 A f + ( } 

where mA, A and A a* are masses and residues of A-isobar and its 

resonance A* respectively. 

Expression ( p4]) is analogous to the contribution from the single-pole term 
appearing in QCD sum rules for correlators in constant external field (see 
I)- 



It is easy to see that the function multiplied by Q 2 /M 2 in (24) changes 
from 1/12 at m 2 — m 2 A —>■ to 1/4 at m\ — m A — > oo. However, taking into 
account the continuum, only contributions from the states with m| — m A ~ 
So — m\ (sq is the continuum threshold) are to be considered. In the region 
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s ^> so, our model of continuum is quite correct, but when s ~ So it is not 
so. Then we see, that nonexponentially suppressed terms will be determined 
by the states with m\ ~ So- Taking m\ — m\ ~ sq — m\ ~ 1.5GeV 2 , 
M 2 ~ lGeV 2 , expression (|24]) can be written in the form 



AaAa * 3 m 2 -ml {1 + + 6) + - } (25) 



where e < 0.1. 



Thus, from (H) and (H) it is seen that when Q 2 /M 2 < 1, and M 



lGeV 2 , the right-hand side of ([19|) is 

\xle-< /M \G M (Q 2 ) + C(Q 2 )M 2 ) e^dz (26) 

The accuracy of (|2~6D is of an order of a few percents (it depends on the 
numerical value of e from (p5|)). It can be shown, that the next terms in 
the expansion (p3|) in powers of Q 2 /M 2 do not change the situation. So, in 
the region Q 2 /M 2 < 1, M 2 ~ lGeV 2 , the right side of the sum rule flT9] ) is 
indeed represented by expression ([26]). 

Let us note that when Q 2 /M 2 > 1, the expression ( PB|) is invalid and we 
cannot separate the contribution of the single-pole terms from the contribu- 
tion of the us double-pole term ~ Gm{Q 2 ) which we are interested in. Thus, 
our sum rule is expected to be valid in the region < Q 2 < lGeV 2 . 

Here we have constructed the right-hand "phenomeno logical" side of the 
sum rule ([BJ) and discussed the region of its applicability. Now, let us pass 
to the calculation of the left, "theoretical" part of our sum rule. 

Using eq.(|7|) and dispersion integral ( |H)D we have obtained the following 
sum rule for Gm- 



r 1/2 e^ 2 dz r dSl r o ds 2 e~^^ pA ^^ Q 

J-l/2 Jo Jo 



2\ 



-1/2 JO JO 7T 2 



+ -a 2 r' 2 dze^ 2 ' M2 (27) 
3 7-1/2 



+ \{2tx) 2 a{i[{ 2 {Q 2 ))M 2 e 



2^-tt /^^2„ 4 ^ 
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|a 2 (l + Q 2 (ffl 4 (Q 2 )))e^ 



-(2vr) 2 ^-n 1 (g 2 )e^ = J(M 2 , s , Q 2 ) [ e^dz 

- m A 1 r+ 1 / 2 q 2 z 2 

= >^le^e- u \-G M {Q 2 ) + C(Q 2 )M 2 / e^cfe 

6 J-l/2 

where e u = +| and Pa(si, s 2 , Q 2 ) is a spectral density: 



^a(*i, s 2 , Q 2 ) = Q 2 ^^ (3(2* + U) + Q 2 (l + 1 k2 ; ) 



where k = (s 1 + s 2 + Q 2 ) - 4siS 2 , U = s 1 + s 2 + Q 2 , a = (2n) 2 \< ^ > 
0.55GeV 3 , A A = (2tt) 2 A a = 2.3GeV 6 , s = AGeV 2 , m A = 1.232GeV, 

in 2 (Q 2 ) = - < ipip > ( -rrr^ — — r — — ) (Gel/) 

Here Q 2 is in units of GeV 2 . 

niW) =(2^(- 9 ^ + m ;) +lD(1 + ^ 







where ^ = 2.3 is residue of the p— meson, m 2 = 0.6GeV 2 is the p— meson 
mass squared, and Sq = 1.5GeV 2 . Numerically in GeV units we have 

n '« 2 ^(-Q^6 +hl < 1 + ^) 

n„(Q 2 ) = in 2 (Q 2 ) 

All these expressions for operators 11^ (Q 2 ) were obtained in [0. In (E7] s = 
4 GeV 2 is continuum threshold, which is numerically equal to the threshold 
of continuum in the sum rules for A-isobar mass. Such a choice for the 
value of the threshold follows from the assumption that there is one and 
the same value of threshold for different structures in the sum rules. Then, 
considering the sum rules for the electric charge (i.e. for the electric form 
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factor Ge{Q 2 ) = 1, Q 2 = 0) it can be shown that they coincide with the 
mass sum rules except for the possibly different values of thresholds. But 
from Ward identity, we know that these sum rules should coincide exactly. 
So the value of the threshold in the sum rules is s Q ~ 4GeV 2 [||. This 
assumption is based on the physical meaning of the quantity s , which as we 
know from the experience is determined mainly by the position of the next 
resonance in corresponding channel. And it would be very surprising if the 
resonance contribution to the electric form factor's channel strongly differs 
from the contribution to the magnetic form factor's channel. 
In the limit So — > oo 

j poo roo s 1 +s 2 . si-s 2 

/ dsi / ds2p A (s 1 ,s 2 ,Q 2 )e~^^ +z ^ T ' (28) 
Jo Jo 



7T 2 

= e u M e (2^(^(1 - Az 2 )) - 3E 3 (fracQ 2 4M 2 (l - Az 2 )) + E 5 (^(l - ^ 

Here E n (z) = $™ x~ n e~ zx dx. 

Now, let us discuss the continuum contribution into the second term of 
the left side of (|27|). The double discontinuity of the corresponding diagram 
could not be calculated due to the reasons discussed in the Section 2 of this 
paper. We do not know the exact formula, analogous to the first term in 
(j2~7D ) but can write an approximate expression, which becomes exact in the 
limit when M 2 — > 00. Moreover, if one substitutes the first term in the left 
side of (p7D by the analogous term 



e„M\l - e-'"-"(l + ^ + Jj))e« a /« 2 x (29) 

C dz[2E ^ {i - 4z2)) - 3£3( ^ (1 - 4z2,) + E ^ {1 - iz2))] 

then one sees that their differences is less then 10% in the region Q 2 < 
lGeV 2 . So, since we shall work in such a region of M 2 , where the continuum 
contribution is less then 20 — 30%, then the uncertainties due to the second 
term in ( p7| ) will be less then a few percents. 

To find Gm{Q 2 ) we should study the following formula to kill nonsup- 
pressed contribution of the nondiagonal transitions 

G M {Q 2 ) = 3e u (l - M 2 ^)A- 2 e m A/^ 2 / (M 2 , S , Q 2 ) 
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at any fixed Q 2 . 

The results obtained for the form factor are depicted in Fig.l. Notice that 
contribution of nondiagonal transitions is very small (about few percents) in 
the sum rule. In the region < Q 2 < 0.8GeV 2 the result obtained may be 
fitted by the following relation 



where p 2 = 0.7 GeV 2 . 

The additive quark model prediction is Ga/(0) = 7 'A^—. 

The analogous sum rule could be also written in the case of f2~-hyperon in 
the limit of SU(3) symmetry after evident interchange e u e s . An accuracy 
of the results obtained is 10 — 20%. 

4 Conclusions 

In the present paper the QCD sum rules for the polarization operator in 
a variable external field are used to calculate the A-isobar magnetic form 
factor in the infra-red < Q 2 < lGeV 2 region. Analogous sum rules can be 
used for calculation of other diagonal hadronic form factors in the infra-red 
region. Notice that this method does not work in the case of nondiagonal 
form factors because it is not possible to separate interesting contribution 
into correlator. There is only one way to do it: to sum all (fcx) ra -terms in the 
expression for nonperturbative propagator (§). It was done only in the case 
Q 2 = [HJ using operator product expansion on a light cone and models for 
the photon wave functions. 
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